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The BCS-BEC crossover and phase diagram for asymmetric nuclear superfluid with pairings in 
isospin 7 = and I — 1 channels are investigated at mean field level, by using a density dependent 
nucleon-nucleon potential. Induced by the in-medium nucleon mass and density dependent coupling 
constants, neutron-proton Cooper pairs could be in BEC state at sufficiently low density, but there 
is no chance for the BEC formation of neutron-neutron and proton-proton pairs at any density and 
asymmetry. We calculate the phase diagram in asymmetry-temperature plane for weakly interacting 
nuclear superfluid, and find that including the 1 = 1 channel changes significantly the phase structure 
at low temperature. There appears a new phase with both 1 = and 7 = 1 pairings at low 
temperature and low asymmetry, and the gapless state in any phase with 7 = 1 pairing is washed 
out and all excited nucleons are fully gapped. 
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I. INTRODUCTION 

As it is well-known, due to the condensate of nucleon- 
nucleon (NN) Cooper pairs at sufRciently low temper- 
ature, a nucleon many-body system, such as a large N 
nucleus or bulk nuclear matter in neutron stars, will be 
in superfluid state with many various interesting phe- 
nomena, like the properties of medium-mass N Z 
nuclei produced at the radioactive nuclear beam facil- 
ities [1, 2], the deuteron formation in medium-energy 
heavy ion collisions [3], and the equation of state of 
neutron stars [4, 5]. Considering the spin and isospin 
degrees of freedom, the NN Cooper pairs should have 
rich inner structure and hence different phase diagrams. 

Recently, two research directions in the study of nu- 
clear matter receive more attention. One is the BCS- 
BEC crossover [3, 6-13]. When nuclear density de- 
creases, the weakly correlated NN BCS state at high 
density may go over to the BEC superfluid of NN bound 
state at lower density. Although the BCS and BEC lim- 
its are physically quite different, the change from BCS 
to BEC was found to be smooth [14-16]. For neutron- 
proton (np) pairs in ^Si Di channel, the chemical 
potential changes sign at a critical density which can 
be regarded as a criterion of the formation of BEC and 
finally approaches to a half of the deuteron binding en- 
ergy at low density limit [3, 7]. Recently, the possible 
BCS-BEC crossover of neutron-neutron (nn) pairs in 
^5*0 channel was studied by examining the spatial struc- 
ture of two correlated neutrons [11, 13] and the density 
and spin correlation functions [10]. It was found that 
a di-neutron BEC state can be formed in symmetric 
nuclear matter at very low neutron density. However, 
when the degree of isospin asymmetry is high, there is 
no such a BEC state in the whole density region. 

The other direction is the possible phase transition 
induced by the mismatch between neutron and pro- 
ton Fermi surfaces [17-22]. When the isospin asym- 
metry becomes sufficiently high, namely when the mis- 
match is comparable with the np pairing gap, the pair- 
ing will be suppressed. A phase transition from the 
BCS state to normal nuclear fiuid is expected at a 



critical isospin asymmetry. However, for an asymmet- 
ric system, besides the BCS state, some other super- 
fiuid states are suggested in condensed matter and nu- 
clear matter, such as the Sarma phase [23] or breached 
pairing phase [24] where the superfluid component 
is breached by the normal component in momentum 
space, the Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) 
phase [25] where a Cooper pair has a total momentum 
and the translational symmetry is spontaneously bro- 
ken, the deformed Fermi surface phase [26, 27] where 
the Fermi surfaces of the two species are deformed 
into ellipsoidal shape and cross to each other so that 
Cooper pairs with zero total momentum could form at 
the cross node, and the phase separation (PS) in real 
space [28, 29] where the normal and superfluid compo- 
nents are inhomogeneously mixed. 

Since nucleons could form Cooper pairs in both 
isospin singlet and triplet channels, it is natural to ask 
a question of how the competition between the isospin 
7 = pairing, namely the np pairing, and 7 = 1 pairing, 
namely the nn and pp pairing, affects the BCS-BEC 
crossover and thermodynamics of nuclear superfluid. 
To have an insight into this question, we perform in 
this paper a mean field analysis for the asymmetric nu- 
clear matter with the above two kinds of pairings. The 
mean field approximation is an effective and successful 
treatment at low temperature or for weakly correlated 
systems. At high temperature or for strongly coupled 
systems, pair fluctuations may become significant [15]. 
We will restrict our study on the BCS-BEC crossover 
at zero temperature and the thermodynamics at high 
density. To simplify numerical calculations. We use a 
density dependent contact potential to describe the NN 
interaction. 

The paper is organized as follows. We present the 
mean field formalism for a general asymmetric nuclear 
matter with isospin 7 = and 7=1 pairings in Section 
II. The possibility of BCS-BEC crossover for the two 
channels is investigated in Section HI, and the signifi- 
cant change in the phase structure of weakly interacting 
nuclear superfluid at flnite temperature is discussed in 
Section IV. We summarize in section V. 
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II. MEAN FIELD FORMALISM 

We adopt the density dependent contact interaction 
(DDCI) developed by Gorrido et. al. [30] to model the 
NN potential, due to its simplicity and validity in pair- 
ing problem. For the purpose of getting qualitative 
conclusions which should not be sensitive to the details 
of the interacting dynamics, the DDCI potential is ac- 
ceptable. The potential is of the form 



F(x,x') 



/ p((x + xO/2)) 



Pa 



5(x-x'), (1) 



where v, rj and 7 are three adjustable parameters, 
/9(x) = p„(x)-f /9p(x) is the nuclear density. Taking suit- 
able values of the parameters, one can reproduce [30] 
the pairing gap A(fci?) as a function of the Fermi mo- 
mentum kp = (37r^p/2)^/'^ in the channels L = 0, 1 = 
1, = ±1, 5 = and L = 0, / = 0, 5 = 1, = 0. We 
will choose in the following numerical calculations the 
parameters [30] rj = 0.45,7 = 0.47,1; = -481 MeVfm^ 
in the 1=1 channel and rj = Q,v = —530 McVfm^ in 
the 1 = channel and the energy cutoff Ec = 60 MeV in 
both channels to regularize the integration. With these 
parameters one must use a density-dependent effective 
nuclcon mass m(p) [30], corresponding to the Gogny 
interaction [31], 



rriQ 
m{p) 



l + ^^5][M/, + 2(i?,-i?, 



X Mce" 



cosh Xc sinh Xc 



(2) 



with Xc = fc|,/i^/2, where mg = 939 MeV is the nuclcon 
mass in vacuum, and /ic, Wc, Bc,Hc, Mc are parameters 
by fitting the Gogny force Dl [32, 33], their values are 
listed in Table I. As shown in Fig. 1, the medium effect 
suppresses the nucleon mass. 



TABLE I: Parameters in the effective nucleon mass (2) by 
fitting the Gogny interaction Dl [32, 33]. 



c ^ic [fm] 


Wc[MeV] 


B^MeV] 


i^4MeV] Mc[MeV] 


1 0.7 


-402.4 


-100.0 


-496.2 -23.56 


2 1.2 


-21.30 


-11.77 


37.27 -68.81 



For a uniform nuclear system, the nuclear density p is 
independent of x, and the DDCI potential is simplified 
to the form of 



F/(x-x') =ff/<5(x-x') 
with the effective coupling constant 

gi = vi[l- r]j{p/pQy'], 



(3) 



(4) 



where / = 0, 1 denote the isospin of pairs. 

The uniform nuclear system with the two-body in- 
teracting potential V/(x — x') and chemical potentials 
Hn and pp for neutrons and protons is controlled by the 
Lagrangian density 



C 



-j d3x'Fi(x-x') [n|(x)nt(x')ni(x')nT(x) +p|(x)p|(x')Pi(x')j5T(x) 
-i J d3xVo(x-x') [nj(x)pj(x')-p|(x)n|(x')] [pi{^')h^{^) - h^{^')p^{x 
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-^50 (n|(x)p|(x) -p|(x)nj(x)) (^j5j,(x)n|(x) - n|(x)pt(x)) , 

I 



(5) 



where fia- and p^- are neutron and proton field opera- 
tors with spin a. We introduce the Cooper pair oper- 
ators A„p(x) = -(go/2)(pi(x)n|(x) - n^(x)p|(x)) = 
A;p(x), A„„(x) = -.gi7i|(x)n|(x) = A,*„(x) and 
App(x) = -5ipx(x)p|(x) = A*p{x.), correspond- 
ing to the channels L = 0,1 = 0,S = l,Sz = 
and L = 0,/ = l,/z^= ±,5'^= 0, and 
keep only their condensates (A„p(x)), (A„„(x)) and 
(App(x)) in the Lagrangian in mean field approxima- 



tion. To incorporate the FFLO state into our study, 
we assume the following forms of the condensates, 
(A„p(x)) = A„pe*(i"+i^)-, (A„„(x)) = A„„e2'i"- 

and (App(x)) = Appe^"*" '', where A„p, A„„ and App 
are constants and can be assumed to be real numbers, 
and q„ and qp are the pair momenta. Obviously, the 
translational symmetry and rotational symmetry in the 
FFLO state are spontaneously broken. Note that, for 
the sake of simplicity, the FFLO state we considered 
here is its simplest pattern, namely the single plane 
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FIG. 1: The effective nucleon mass as a function of nuclear 
density, calculated from the Gogny force Dl. mo and po are, 
respectively, nucleon mass in vacuum and normal nuclear 
density. 



wave FFLO state or the so-called FF state. 

The partition function Z which is the key quantity 
of a thermodynamic system can be calculated by path 
integral, 



[dn„] [dp„] [dnl] [dp^] exp 




d^xC 



(6) 

By performing a gauge transformation for the nucleon 



fields fia 



and Pa = e ^'^'''^Pcr which keeps 



the functional measure in Eq.(6) invariant, the path in- 
tegral over fig- and p^ can be easily done and we obtain 
the mean field thermodynamic potential 
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90 



■ A2 
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d^k 
(27r)3 



TrlnG-l(^cJ„k), (7) 



where lOi, = {21/ -\- 1)'kT with i/ G Z is the fermion 
frequency, and G is the Nambu-Gorkov propagator, 
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\ ^nn ^np 

with the definition = (k±qi)^/(2m) — /i^ for i = n,p. 
It is worthy noting that, if the relative momentum 
~ ^p is large enough, the uniform superfluid would 
be unstable due to the stratification of the superfluid 
components characterized by A„„ and App, in analo- 
gous to the multi-component BEG in condensed matter 
physics. To avoid such a dynamic instability [34], we 
choose Qn = Qp = q. After computing the frequency 
summation and trace in (7), the thermodynamic poten- 



tial can be expressed in terms of quasi-particles, 

d^k 



2A2 A2 + 

^ np nn ' pp 

90 91 

'si 



(2n) 



(9) 



where are the quasi-particle energies 



with e±, 6e and ea defined as 



2e± — 



A^ + A^ 

nn ' np 



± (e+Ep 



T5e 



pp 



(10) 



A 



np } 



5e^{e+-en)/2^{e+-,-)/2, (11) 
4 = A2p[(e+ - e+)(e- - ) + (App - A„„)2]. 

From the thermodynamic potential, we derive the 
neutron and proton number densities. 



Pn 



dy.n 



"^■^ d^p 



or 



9/i ' 



6p 



(13) 



where p, = {pn + Pp)/2 is the averaged nucleon chemical 
potential, Sfi ~ (fin — Pp)/2 the mismatch between ^„ 
and /ip, and 6p = p„ — pp the associated difference 
in number densities. We will focus on neutron rich 
nuclear matter with 6p > and Sp > 0, since it is 
coincident with the environment of neutron stars and 
some heavy nuclei. For the sake of convenience, we 
define the relative density asymmetry a = Sp/p. 

The condensates and the FFLO momentum are de- 
termined by the gap equations. 



on 

dA 



= 0, 



np 



dAnr. 



= 0, 



dn dn 

= 0, - = 0. (14) 



OA. 



pp 



With the above coupled number equations (13) and 
gap equations (14). we can solve the condensates 
A„p, Ann, App, FFLO momentum q and chemical po- 
tentials fi, Sp as functions of temperature T and den- 
sities p, Sp. The ground state of the system is speci- 
fied by the solution corresponding to the global mini- 
mum of the free energy which is related to the ther- 
modynamic potential by a Legendre transformation, 

= n + PnPn + PpPp = n + I^Lp + SpSp. 



III. BCS-BEC CROSSOVER AT ZERO 
TEMPERATURE 

The mean field approximation has been successfully 
used to study the BGS-BEG crossover in nuclear matter 
with isospin / = [3, 6-9] or / = 1 [10-13] pairings. 
We extend in this section the study to a more general 
case with both / = and 7 = 1 pairings. To guarantee 
the validity of mean field treatment, we focus on the 
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supcrfluid at zero temperature. Since the FFLO state 
is unstable in strong coupling region [35], we will not 
take into account the FFLO state in the study of BCS- 
BEC crossover. The temperature effect and the FFLO 
state will be considered in Section IV. 

To examine how the NN correlation changes with 
nuclear density, we define the normalized Cooper pair 
wave function as 

V'y(r) = C{BCS\al^{^)a]^{^ + r)\BCS) 

with i,j = n,p, where a]^ is the nucleon creation oper- 
ator, and ^pij(]s.) the wave function in momentum space 
or the anomalous density 

V'y(k) = {BCS\al^ik)ali-k)\BCS). (16) 

The density distribution functions ni(k) is defined 
through 

'^^W = I 5](i3C5|al(k)a.,(k)|i?C^). (17) 

(T 

In order to describe the BCS-BEC crossover quan- 
titatively, it is convenient to introduce the following 
characteristic quantitics[ll, 13]: 

1) The distribution function r^|V'ij(r)p which is the 
probability to find a pair of nucleons i and j with dis- 
tance r in between. For the BCS-BEC crossover in- 
duced by changing coupling constant at fixed density 
such as in cold atom gas, the probability would dis- 
tribute in a wide space region in the weakly coupled 
BCS state but peak sharply at a small r in the strongly 
coupled BEC state. However, for nuclear superfluid, 
the BCS-BEC crossover is induced by changing the nu- 
clear density. The density dependence of any physical 
quantity of the system is reflected in two aspects: the 
explicit or direct p dependence through the coupling be- 
tween the gap and number equations and the indirect 
p dependence through the effective coupling constant 
gi{p) and in- medium nucleon mass m{p). When the 
density drops down, the indirect p dependence is the 
driving force for the crossover from BCS to BEC, but 
the explicit p dependence makes the system dilute and 
then slows down or even rejects the crossover. The 
balance between these two opposite aspects controls 
the BCS-BEC crossover and will probably change the 
monotonous shrinking behavior of the probability. 

2) The root-mean-square radius of the Cooper pair 
= \/ with {r'^)ij = J (fr r'^\ij},j{r)\'^ which 

characterizes the size of the Cooper pair. When the 
density is fixed, ^ij is expected to be large in the weakly 
coupled BCS region and small in the strongly coupled 
BEC region. However, for nuclear superfluid with de- 
creasing nuclear density from BCS to BEC, itself 
can no longer describe the BCS-BEC crossover. In this 
case, the BCS and BEC states are defined in the sense 
of overlapping degree of pair wave functions. At high 
density the size of a pair can be small but the pairs 
may overlap strongly, while at low density the size of 
a pair is probably large but the overlapping of pairs 



becomes weak. To represent the density dependence of 
the overlapping, we can compare with the averaged 

— 1/3 

distance c?y = between any two nucleons with 

Pnp = p/2, Pnn = Pn and = pp. For $y > 
the overlapping is strong, and the ij— pair can be in- 
terpreted as an extended BCS pair, while for <C dij 
the overlapping becomes weak, and the pair should be 
considered as a compact BEC pair, i.e., a boson- like 
bound state. 

3) The probability P(dy ) = An J^'' dr r'^\ij,j{r)\^ of 
finding an ij— Cooper pair within the averaged distance 
dij between the nucleons i and j. Note that the proba- 
bility is density dependent and can be used to describe 
the overlapping degree of pairs. It should be very close 
to 1 in the strongly coupled BEC state and clearly less 
than 1 when the superfluid is in the BCS state. 

4) The s-wave scattering length Uij which relates the 
coupling constant gj to the low energy limit of the T- 
matrix for nucleons i and j in vacuum. 



with ek = k2/(2m) and 5„p = go,5n« = 9pp = 9i- 
In the BCS region is negative, representing the at- 
tractive interaction between nucleons. However, in the 
BEC region it should be positive to preserve the stabil- 
ity of the two-body bound state. Therefore, the change 
of ttij from negative to positive value can be considered 
as a signature of the BCS-BEC crossover. 

5) The scaled condensate /e^ with the Fermi energy 
defined as 4 = {^i'f I i'^rn) , where kj = kp, kj" = 
knF and k^ = kpp are the corresponding Fermi mo- 
menta. In the case with only / = pairs, the mean field 
gap equation and number equation satisfy the so-called 
universality: The scaled quantities like A/ep and p/f-F 
are only functions of the effective coupling l/{kFa) [14]. 
In the strongly coupled BEC region, one should expect 
a large A/ef, but in the weakly interacting BCS state 
it will be less than 1 . While the universality is explicitly 
broken in our general case with both / = and I ~ 1 
pairings, we still take the scaled quantities to describe 
the BCS-BEC crossover. 

6) The effective chemical potential pij/sp with pnp = 
fJ-, P-nn = fJ-n and ppp = Pp. Due to the definition, 
Pij / Ep is exactly equal to 1 in the BCS limit and the 
order of 1 in the weakly coupled BCS region. In the case 
of only one kind of pairs, like symmetric nuclear matter 
and neutron matter discussed below, p/ef is negative 
in the BEC region, the absolute value \p/ef\ becomes 
larger and larger as the system goes deeper and deeper 
into the BEC region, and 2p can be viewed as the bind- 
ing energy of the bound state in the BEC limit when 
the system approaches to vanishing density. That is the 
reason why the change of the sign of p is normally con- 
sidered as an unambiguous criterion of the formation of 
BEC, at least at mean field level. However, in general 
asymmetric nuclear matter with np, nn and pp pairs, 
the correlations among different pairs will qualitatively 
change this conclusion, see the detailed discussion in 
the following. 
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A. Symmetric Nuclear Matter: q = 



Let us first consider the symmetric nuclear matter 

The number and 



with 6p = 0, 6fi = and A„„ 
gap equations in this case arc largely simplified as 



Pn 

ri(k) 



Pp 



(2^ 



n(k), 



n„(k) = np(k) = 2 ' ^ 



d^k 1 



d^k 1 



(27r)3 2ii;k' 



(19) 



where ek = k'^/(2TO) 
lation and i?k = 



/i is the particle dispersion re- 
A^j„ the quasi-particle 



A2 

np 



trivial solutions A„p 7^ 0, A„„ = 
0, A„„ = App 7^ 0. Considering 



dispersion. The gap equations have two types of non- 

= and A„p = 
gi\ at any den- 
sity, only the solution with A„p 7^ corresponds to the 
ground state. This can also be verified by comparing 
the free energies for the two solutions. Since there is 
no experimental evidence for np pairing at finite nu- 
clear density, it is necessary to note that the condition 
Iffol > |.9i| is in principle an assumption in our treat- 
ment. From a more self-consistent Greens function ap- 
proach [36] , the tensor correlations in the / = channel 
may yield a kind of pair correlation which is different 
from the one observed in solving gap equations. 

To simplify the notification, we write in short ep = 
Cp and d = dnp in the following. Solving the coupled 
number and gap equations for the np channel, we ob- 
tain the scaled gap A„p{p)/eF and chemical potential 
fi{p)/eF as functions of nuclear density. 

Defining the anomalous density '0„p(k) = A„p/(2i?k) 
and substituting it into the number and gap equations 
give the Schrodinger-like equation 



k2 

-0np(k) 



(1 - 2nk).9o / T^-TT^npC*') = 2^V'np(k) 

(20) 

for the anomalous density. In the limit of vanish- 
ing density, rik 0, the equation goes over into the 
Schrodinger equation for the np bound state, and the 
chemical potential 2^ plays the role of binding energy. 
Since such a bound state should be a boson, one expects 
that, at sufficiently low density and low temperature 
the symmetric nuclear matter is in the BEC phase. 

Making the Fourier transformation of the anomalous 
density '0„p(k) which is determined by solving the gap 
and number equations, we obtain the wave function 
i^npi'^) in coordinate space. The behavior of the prob- 
ability distribution r^|'!/)„p(r)p as a function of the rel- 
ative distance r between the pair partners is shown in 
Fig. 2. The spatial extension and profile of the prob- 
ability depend strongly on the density. Near the nor- 
mal density, the probability function is spatially ex- 
tended and behaves like the well-known BCS expres- 
sion [37] ij^npir) ~ KQ{r/TT£_np)s.va.{kFr)/{kFr). The os- 
cillation induced by the sharp Fermi surface at high 
density is well-known and called Friedel oscillation [38] 




FIG. 2: The probability r^|?/'„p(r)l^ as a function of the 
relative distance r between the paired neutron and proton 
at different nuclear density p in symmetric nuclear matter. 



which is widely discussed in nuclear matter and recently 
extended to quark matter [39]. The Fermi surface is 
strictly defined only in normal state without any con- 
densate. In the BCS region at high density, the con- 
densate is small and there is still an approximate Fermi 
surface. However, in the BEC region at low density, the 
Fermi surface is already very low and further deformed 
by the large condensate. The number density as a func- 
tion of momentum is plotted in Fig. 3. It is clear that 
the Fermi surface is well-defined only at high density. 
This is the reason why the Friedel oscillation almost 
disappears at low density. 
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FIG. 3: The nuclear density distribution n(fc) as a func- 
tion of momentum at fixed nuclear density p in symmetric 
nuclear matter. fc_F is the Fermi momentum. 



It is necessary to note that the probability 
''^|'0np('')|^ does not shrink monotonously with decreas- 
ing number density, see Fig. 2. For a very dilute system, 
the two partners of a pair can have a large distance be- 
tween them, the pair wave function will be spatially 
extended, and the probability can then distribute in a 
wide space region. While the strong Friedel oscillation 
at high density and its weakening at low density can be 
considered to distinguish the BCS state from the BEC 
state, the space extension of the probability itself can 
not be used to characterize the BCS-BEC crossover in 
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nuclear supcrfluid. 

To clearly describe the BCS-BEC crossover induced 
by the change in nuclear density, we now calculate the 
scaled root-mean-square radius ^„p/d, the scaled chem- 
ical potential ^/cf, the scaled gap parameter Anp/epi 
the effective scattering length l/{kpanp) and the prob- 
ability Pnp{d). inp/d and ii/ep are presented in Fig. 4 
as functions of nuclear density. While ^„p itself is not 
a monotonous function of p, the scaled one goes up 
monotonously with increasing density. The right verti- 
cal line at p/po ~ 0.4 indicates the position of S,np/d — 1 
which can be used to separate the BCS region with 
sharply increasing S,np/d from the region with slightly 
changing and small S,np/d. The scaled chemical po- 
tential is 1 in the BCS limit, then drops down with 
decreasing density, and becomes negative at very small 
density. The position where p/cf approaches to zero 
is indicated by the left vertical line at p/po ^ 0.003 
which is, from the definition, considered to distinguish 
the BEC region with negative p from the other region 
with positive p,. Therefore, considering both £,np/d and 
p/ep, the BCS and BEC states are, respectively, lo- 
cated at p/ po > 0.4 and p/ po < 0.003, and the crossover 
from BCS to BEC is in between the two vertical lines. 
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FIG. 4: The scaled root-mean-square radius (,np/d and 
scaled chemical potential p/ep as functions of nuclear den- 
sity in symmetric nuclear matter, ef is the nucleon Fermi 
energy and d the averaged distance between two nucleons. 
The right and left vertical lines which separate the BCS, 
BEC and crossover regions are, respectively, determined by 
the conditions ^„p/d — 1 and p/ep = 0. 



Pnp{d) decreases monotonously with increasing density. 
It approaches to 1 at low density which means strong 
correlation of pairs in the BEC state and becomes very 
small at high density which indicates a weak correlation 
of pairs in the BCS state. The two vertical straight lines 
in Fig. 5 and Fig. 6 which separate the BCS, BEC and 
crossover regions are still characterized by the scaled 
root-mean-squarc radius and scaled chemical potential. 
By comparing the above three figures, we can see that 
Cnp/d, p/f-F, Anp/ep, l/{kpanp) and P„p((i) can self- 
consistently describe the BCS-BEC crossover in sym- 
metric nuclear supcrfluid. 
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FIG. 5; The scaled condensate Anp/ep and effective scat- 
tering length l/{kpanp) as functions of nuclear density in 
symmetric nuclear matter. a„p is the neutron-proton scat- 
tering length. 
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In Fig. 5 wc show the scaled condensate Anp/ep and 
the effective s-wave scattering length 1/ (kpanp). Both 
are monotonous functions of p, while again the gap pa- 
rameter A„p itself does not behave monotonously. The 
scaled gap is large at low density and small at high 
density and approaches to zero in the BCS limit. The 
effective scattering length drops down with increasing 
density and becomes negative at high density. The den- 
sity dependence of A„p/eF and \/ {kpanp) agrees well 
with the normal understanding of both BCS state with 
small gap and negative scattering length and BEC state 
with large gap and large and positive scattering length. 

The probability Pnp{d) of finding a pair with relative 
distance r < d between the paired neutron and pro- 
ton is plotted in Fig. 6. Again, while the probability 
Pupif) is not a monotonous function of p for fixed r, 



FIG. 6: The probability Pnp{d) of finding a np-pair with 
relative distance r < d between the paired neutron and 
proton as a function of nuclear density in symmetric nuclear 
matter. 



B. Neutron Matter: a = 1 

Neutron matter, as fully asymmetric nuclear mat- 
ter, is important for physics of neutron stars. In 
this case, there are a = l,p„ = p, Pp — and 
Anp = App = 0. By solving the coupled number equa- 
tion and gap equation for p„ and A„„ which are ex- 
actly the same as Eq.(19) when we replace ek and £'k 
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k A2„, we 



by Enk = k2/(2m) - ^„ and £;„k 
obtain the neutron chemical potential /i„ and nn pair 
condensate A„„ as functions of nuclear density. The 
anomalous density ?/;„„(k) = A„„/(2£'„k) satisfies the 
similar Schrodinger-like equation (20), and in the low 
density limit 2/i„ plays the role of binding energy of the 
possible di-neutron bound state. 

In comparison with the symmetric nuclear matter, 
the probabihty r^|V'nn(r)P and the density distribution 
n„(k) at different total number density are very simi- 
lar to r^|?/'„p(r)p and ri(k) shown in Fig. 2 and Fig. 3, 
but the scaled root-mean-square radius ^nn/dnm chem- 
ical potential Hn/^nF, gap parameter A„„/e„i? and the 
effective scattering length l/(fc„j?a„„) behave very dif- 
ferently in the low density region. In Fig. 7 and Fig. 8, 
the vertical straight line located at p = O.lSpo is de- 
termined by the condition ^nn/dnn ~ 1, which indi- 
cates the BCS region at higher density. However, one 
can not define the BEC region through the definition 
/i„ = 0, since /i„ is always positive in the whole den- 
sity region, and correspondingly the scaled condensate 
A„„/e„F is still small and the effective scattering length 
l/(fcnFinn) is Still positivc at extremely low density. 
All of these characteristics indicate that no boson de- 
gree of freedom emerges in neutron matter at low den- 
sity. 
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FIG. 7: The scaled root-mean-square radius S^nn/dnn and 
neutron chemical potential [inj^nF as functions of nuclear 
density in neutron matter. The vertical line represents the 
BCS boundary. ci„„ is the averaged distance between two 
neutrons and e„F the neutron Fermi energy. 



One can understand the reason why there is no BEC 
in neutron matter in terms of the density dependent 
coupling constant gi and the Fermi surface knP ■ For the 
DDCI potential parameters chosen by fitting the pair- 
ing gap versus Fermi momentum, the effective coupling 
in / = I channel is weaker than the one in / = channel 
at any density, |(7i(p)| < |3o(p)|i and on the other hand, 
the neutron density p„ = p in neutron matter is two 
times the neutron or proton density pn = Pp = p/2 m 
symmetric matter. From this comparison, while BEC 
can form in symmetric nuclear matter, its formation in 
a dense neutron matter with weak coupling becomes 
difficult and even impossible. This explains also why 
the BCS boundary shifts from p/po = 0.4 in symmetric 
matter to 0.15 in neutron matter. 
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FIG. 8: The scaled condensate A„„/eF and effective scat- 
tering length l/(kF functions of nuclear density 
in neutron matter. a„„ is the neutron-neutron scattering 
length. 



C. Asymmetric Nuclear Matter: < a < 1 

While only np condensate in symmetric nuclear mat- 
ter with a ~ Q and nn condensate in neutron mat- 
ter with a = 1 can survive, there may exist np, nn 
and pp condensates in asymmetric nuclear matter with 
< a < 1. In this general case, the three condensates 
A„p, Ann and App and the chemical potentials p and 
Sp as functions of nuclear density p and asymmetry Sp 
are calculated by solving the coupled three gap equa- 
tions and two number equations, and the ground state 
of the system is determined by comparing the corre- 
sponding free energies. 



0.04 




FIG. 9: The probabilities r^\il'np{r)\'^ , r^\ipnnir)f and 
r^|i/)pp(r)p as functions of the relative distance r between 
the pair partners at different nuclear density in nuclear mat- 
ter with asymmetry a = 0.1. 



The three probabilities r'^\ipnp{i")\^ , r^|i/)„„(r)p and 
r^|'!/)pp(r)p for np, nn and Cooper pairs are shown in 
Fig. 9 at different nuclear density and fixed asymmetry 
a ~ 0.1. At high density p/po = 0.6, all the proba- 
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bilities are spatially extended with strong Friedel oscil- 
lations, indicating that all paired nucleons are weakly 
correlated and the ground state is in the BCS phase 
with all three condensates. Similar to the case in sym- 
metric nuclear matter, with decreasing nuclear density 
p the np pairing probability shrinks first but slightly 
expands again at extremely low density. The very sur- 
prising feature is that the nn and pp pairing proba- 
bility are spatially expanded in a wide region at low 
density and even approximately r-independent at ex- 
tremely low density. This is a strong hint that there 
are no BEC states of nn and pp pairings in general 
asymmetric nuclear matter. 
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FIG. 10: The scaled root-mean-square radius £,ij/dij and 
scaled chemical potential in/ep as functions of nuclear den- 
sity in nuclear matter with asymmetry a = 0.1. 



To check if the system could reach BEC state at 
low density, we calculate the scaled root-mean-square 
radius £,ij/dij and chemical potential fii/ep and show 
them in Fig. 10 as functions of nuclear density at fixed 
asymmetry a = 0.1. The BCS region for the ij pairing 
is defined through the condition ^ij/dij = 1. The BCS 
boundary is roughly at p/po = 0.35 for np pairs, 0.14 for 
pp pairs and 0.17 for nn pairs. Below the BCS bound- 
ary the system is regarded as a strongly correlated su- 
perfluid. As mentioned above, for a system with only 
one kind of pairings such as symmetric nuclear matter, 
the strongly correlated system will go into the BEC 
state when the chemical potential becomes negative. 
With this judgement, the np Cooper pairs could reach 
BEC state at the critical density p/po = 0.002 where 
the chemical potential n changes sign. For the nn pairs, 
they could never form di-neutron bound state, since the 
neutron chemical potential is always positive, although 
it drops down as density decreases. The interesting 
phenomenon is for the pp pairing. At p/po = 0.03, 
the proton chemical potential changes sign, even earlier 
than the change for the np pairing! Does the negative 
Pp here mean the formation of pp BEC state? Is it 
inconsistent with the flat structure of the probability 
shown in Fig. 9? To answer this questions, we should 
note that for a general asymmetric system, the three 
condensates are strongly coupled and for any anoma- 
lous density V'u(k) = Ay/(2£'ik) there is no simple 
Schrodinger-like equation Eq.(20), and the correspond- 
ing chemical potential pi at low density limit can not 



be identified as a half of the binding energy of the di- 
nucleon bound state. Therefore, there is no longer a 
definite relation between the sign change of pi and the 
BEC formation. 

To verify whether the BEC state is reached, we cal- 
culate the effective s-wave scattering length l/{kpaij), 
the scaled condensate Aij/eip, and the probability 
Pij (d) as functions of nuclear density at fixed asym- 
metry a ~ 0.1. The results are presented in Fig. 11 
and Fig. 12. The scattering length in / = channel 
changes its sign at about p/po — 0.07, but the other 
two lengthes in / = 1 channel remain negative in the 
whole density region. Correspondingly, the scaled con- 
densate Anp/ep becomes much larger than 1 at low 
density but A„„/e„i? and App/epp are still very small 
even at extremely low density, and Pnp{d) is almost 
equal to 1 but Pnn{d) and Ppp{d) are clearly less than 
1 at low density. Therefore, we can safely say that, the 
np BEC state is reached at low density, but there are 
no nn and pp BEC states in general asymmetric nuclear 
matter. 
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FIG. 11: The scaled condensates Anp/ep, Ann/^nF, 
App/epF and scattering lengthes l/{kFa„p), l/ikFCinn), 
l/ikpapp) as functions of nuclear density in nuclear mat- 
ter with asymmetry a = 0.1. 
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FIG. 12: The probability Pij{d) as a function of nuclear 
density in nuclear matter with asymmetry a = 0.1. 



The behavior of the nuclear superfluid in general 
case with < a < 1 is controlled by the two con- 
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straints \gi{p)\ < |.9o(/o)| and pp = (1 - a)/2p < p„ = 
(1 + a)/2p < p. For np pairings, the effective coupling 
constant go is density independent in our case, and the 
crossover from BCS to BEC is characterized by the den- 
sity dependent nucleon mass (2). For nn or pp pairings, 
the increase of the density leads to a high Fermi surface 
and then makes the pairing easy, but on the other hand 
the coupling constant drops down with increasing den- 
sity. It is the balance between these two effects that the 
probability P„„((i) or Ppp[d) is not a monotonous func- 
tion of density but there exists a valley structure at low 
density. Since the coupling constant gi is the same for 
nn and pp pairings but the Fermi surface for nn pair- 
ing is higher than that for pp pairing, pp pairing is most 
impossible to form a BEC state, that is the reason why 
the scaled condensate and probability for pp pairing 
are less than the corresponding values for nn pairing, 
see Fig. 11 and Fig. 12. When we consider Coulomb ef- 
fect which is neglected in our treatment and which will 
break the relation = gpp, the much stronger screen- 
ing effect for pp interaction leads to \gpp\ < \gnn\- As a 
result, pp pairing will become further impossible to be 
in BEC state. 



D. Phase Diagram in p — a Plane 

The phase diagram of asymmetric nuclear matter in 
the p ^ a plane at zero temperature is presented in 
Fig. 13. For a = only A„p could exist, and for a — 1 
only Ann survives. These two limits have been dis- 
cussed in III A and III B, and we will in this sub-section 
consider the phase structure in the region of < a < 1. 

Due to the mismatch between neutron and proton 
Fermi surfaces induced by the asymmetry a 0, there 
should be no np pairing when the mismatch is large 
enough. Since the coupling constant in / = channel is 
density independent and the pairing is dynamically con- 
trolled only by the in- medium mass, the critical asym- 
metry anp{p) for the BCS superfluid increases with 
decreasing density. When the system enters the su- 
perfluid state, the correlation becomes more and more 
strong with decreasing density and the system will go 
into the crossover region at a critical density defined by 
^np/d = 1. When the density further decreases, the np 
pairing starts to be in the BEC state at another critical 
density defined hy p = 0. Corresponding to the above 
two critical densities, the two boundaries of the BCS- 
BEC crossover in Fig. 13 are indicated respectively by 
"«p(p) and alpip). At higher asymmetry, a^p and alp 
coincide and the nuclear matter goes directly from the 
normal state into the strongly correlated superfluid. 

Different from the / = channel, the pairing in / = 1 
channel has no Fermi surface mismatch and the super- 
fluid can survive at any asymmetry when the density 
is not high enough. In Fig. 13 the I — 1 superfluid ex- 
ists in the whole p — a plane. As we mentioned above, 
at q; = 0.1 the critical density for nn pairing to go into 
the strongly coupled region is larger than the one for pp 
pairing. This is true for any asymmetry. In Fig. 13 the 
boundary a^nip) defined by inn/dn = 1 for nn pairing 
is always on the right-hand side of the one app{p) de- 
fined by S,pp/dp = 1 for pp pairing. It is interesting that 



the critical density for nn pairing is almost a constant 
for any asymmetry. There is no room for the BEC state 
of nn or pp pairing at any asymmetry < a < 1. 

Note that there exist three jumps in the phase dia- 
gram. The jump from normal state at p = to the np 
pairing BEC state at p 7^ 0, the jump from np pairing 
at a = to the np, nn and pp pairings at a 7^ 0, and 
the jump from nn pairing at a = 1 to np, nn and pp 
pairings at a < 1. 




iogio(p/po) 

FIG. 13: The phase diagram for asymmetric nuclear mat- 
ter in p — a plane at zero temperature, a^p is the critical 
asymmetry for np pairing BCS superfluid, a^p and a^p are 
the two boundaries of np pairing BCS-BEC crossover, and 
o-nn {o.pp) is the boundary for nn (pp) pairing to enter the 
strongly coupled region. 



IV. THERMODYNAMICS AT HIGH DENSITY 

The Cooper pairs could be broken by both large mis- 
match between the Fermi surfaces of the partner nu- 
cleons and large thermal fluctuation. The former is 
characterized by the asymmetry a, and the latter is 
described by temperature T. We in this section exam- 
ine the thermal and asymmetric effect on the nuclear 
matter in weak coupling region, i.e., in the BCS region 
where the mean field analysis is still reliable. From the 
above calculation at zero temperature, this region is 
around and below the normal nuclear density pg. Our 
purpose is to determine the phase diagram in a — T 
plane at a fixed density. 

We must emphasize that the asymmetry may lead to 
the so-called phase separating instability, which means 
that the ground state may favor the spatial separation 
of different phases. To study such an instability we 
define the nuclear density susceptibility matrix x with 
elements Xij = dpi/dpj with i,j = n,p. For an uni- 
form matter, x is always positively defined, and x < 
can be used as an indication of the appearance of phase 
separation (One should note that x < is only a nec- 
essary but not sufficient condition for the occurrence of 
phase separation) . In the following discussion the term 
phase separation is refer to the state with x < 0. In a 
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practical manner, x can be computed through [40] 

Xij - 
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FIG. 14: The phase diagram of a mismatched nuclear super- 
fluid at fixed nuclear densityp/po ~ 0.6 in a — T plane. The 
upper panel is for the familiar superfluid with only 7 = 
pairing, and in the lower panel the 7 = 1 pairing is included 
as well. The shadowed regions indicate the gapless super- 
fiuid. Ao = 7.7 MeV is the condensate in 7 = channel at 
a = and T = 0. 



By calculating the coupled gap and number equa- 
tions at finite temperature, we obtain all the possible 
homogeneous phases and inhomogeneous FFLO phase. 
From the comparison of their free energies, we then 
extract the lowest one at fixed p, a and T. Finally 
we investigate the stability of the system against the 
number fluctuations by computing the number suscep- 
tibility matrix x, the state with non positive-definite 



X may be a phase separation. In Fig. 14 we show the 
phase diagram in a — T plane at fixed nuclear density 
p/po = 0.6. 

The phase diagram without 7=1 pairing is pre- 
sented in the upper panel. The system is in nor- 
mal state when the asymmetry or temperature is high 
enough, and the superfluid is in homogeneous state 
at high temperature and FFLO state at low tempera- 
ture. However, the number susceptibility in the region 
of low temperature and low number asymmetry is not 
positive-definite, the FFLO state in this region is there- 
fore unstable against the number fluctuations, and the 
ground state is probably an inhomogeneous mixture of 
the BCS superfluid and normal nuclear fluid. The shad- 
owed region is the gapless superfluid with Sp > A„p 
where the energy gap to excite quasi-nucleons is zero 
and the system may be sensitive to the thermal and 
quantum fluctuations. 

The phase diagram with both 7 = and 7=1 pair- 
ings is shown in the lower panel of Fig. 14. Besides the 
familiar phase with only 7 = pairing (A„j, 7^ 0, A„„ = 
App = 0) and the expected phases with only 7=1 pair- 
ing (A„„, App 7^ 0, A = and A„„ 7^ 0, A„p = App = 
0), there appears a new phase where the two kinds of 
pairings coexist (A„p, A„„, App 7^ 0). In this new phase 
the FFLO momentum is zero and the number suscep- 
tibility is negative, x < 0- Therefore, the homoge- 
neous superfluid in this region is unstable against num- 
ber fluctuations, and the ground state is probably an 
inhomogeneous mixture of these three superfluid com- 
ponents. In the familiar phase with only 7 = pairing, 
there remains a stable FFLO region and an unstable 
FFLO triangle where the number susceptibility is neg- 
ative and the system may be in the phase separation of 
the BCS superfluid and normal nuclear fluid. The gap- 
less state appears only in the 7 = pairing superfluid, 
and in the region with 7 = 1 pairing all the nucleons 
arc fully gapped. Since the Fermi surface for nn pair- 
ing is higher than the one for pp pairing, the critical 
temperature to melt the nn condensate is higher than 
that for melting the pp condensate, and the difference 
between the two increases with increasing asymmetry. 
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V. SUMMARY 

We have investigated the phase structure of isospin 
asymmetric nuclear superfluid with pairings in both 
/ = and 1=1 channels in the frame of the den- 
sity dependent contact potential. We calculated at 
zero temperature and in mean field approximation the 
pair wave functions, pair condensates, nuclcon chemical 
potentials and effective couplings which are normally 
considered as characteristic quantities describing BCS- 
BEC crossover, and found that in general asymmetric 
nuclear matter only the np pair could form true bound 
state at extremely low density, and nn and pp pairs, on 
the other hand, could never form bound state at any 
density and asymmetry. We also studied the phase di- 
agram for weakly coupled nuclear superfluid at finite 
temperature. Since the attractive interaction for 1=1 



pairing is weaker than the one for 7 = pairing, the 
inclusion of nn and pp pairings changes significantly 
the conventional phase diagram with only np pairing 
only at low temperature. For systems with fixed nu- 
clear density, the two kinds of pairings can coexist at 
low temperature and low number asymmetry. By cal- 
culating the number susceptibility, this new phase is 
not in the FFLO state but probably an inhomogcneous 
mixture of the np, nn and pp superfluid components. 
In any region with 1 = 1 pairing, the interesting gap- 
less superfluid is washed out and all quasi- nucleons are 
fully gapped. 
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